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I. INTRODUCTION 

Over the last decade, a lot of progress has been made 
on the solution of toy models belonging to the Gross- 
Neveu family [l|. They describe N species of massive or 
massless fcrmions in 1 + 1 dimensions interacting via a 
scalar or pscudoscalar four-fcrmion interaction: 
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C = V>(i7 M d M -TO)^ + ^(V^r + — (iph V) ■ (1) state 



Summation over flavor indices is implied and we abbre- 
viate ijjip = ip^i 1 ^ ■ For G 2 = and m = we re- 
cover the original Gross-Neveu model (GN) with discrete 
chiral symmetry ip — > j 5 tp pj. For g 2 = G 2 the theory 
possesses a continuous chiral symmetry ip — > e 17 e ip and 
corresponds to a 1 + 1 dimensional version of the Nambu- 
Jona-Lasinio model (xGN) Q. Over the past years, pre- 
vious work on the phase diagrams of those models p|, H| 
has been extended significantly and revised phase dia- 
grams were proposed that respect the particle content of 
the theories. It was found that all of the above models ex- 
hibit crystalline phases, for many of which analytical so- 
lutions could be obtained The GN model features 
a kink-antikink crystal at high density whereas the ther- 
modynamically preferred ground state of the xGN model 
below a transition temperature is a "chiral spiral" — a 
helical chiral condensate whose amplitude is determined 
by the temperature and spatial period by the chemical 
potential [5[. The phase transition is reminiscent of the 
Peierls transition in Id metals with the first band gap 
opening up at the Fermi surface [ll|. While those mod- 
els provide a rich playground for the study of rclativistic 
field theories, they were not thought to have any appli- 
cation to reality due to their lower dimensionality. 

In recent years, however, the expertise on those toy 
models has been used to study problems related to the 
ground state of QCD at finite density. Thereby, 1 + 1 di- 



mensional models arise in the form of effective field theo- 
ries when phases with lower-dimensional modulations of 
the chiral condensate are considered. Such modulations 
can occur due to a strong external field [l2j or can be 
induced by the Fermi surface HH . Examples include 
work by Shuster and Son [l3| who — based on a dimen- 
sional reduction onto a variant of the chirally invariant 
Thirring model — refuted the possibility of the large N c 
Deryagin-Grigoriev-Rubakov (DGR) chiral wave ground 



15| for N c = 3. More recently, Pcicrls-like insta- 
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bilities were discussed in the context of the Quarkyonic 
Phase of QCD [n| or in the presence of a very strong 
magnetic field, where they were named " Quarkyonic Chi- 
ral Spirals" [3, [l7| and "Chiral Magnetic Spirals" [Hj], 
respectively. The exact solution of the massive GN model 
was used by Nickel in order to investigate the phase di- 
agram of the massive Nambu-Jona-Lasinio model when 
restricting to one-dimensional condensates by mapping 
the energy spectrum of the model onto the spectrum of 
the GN model fl8|. Further examples include work by 
Bietenholz et al. (l9j . 

Since the models in question are very involved, it is 
desirable to study examples of dimensional reduction iso- 
lated from the context of QCD. For this purpose we in- 
vestigate dimensional reduction in nonrelativistic quan- 
tum field theories in this paper. We establish the equiv- 
alence of BCS theory of superconductivity to the mass- 
less Gross-Neveu model with discrete chiral symmetry 
on a mean-field level, provided that the chiral conden- 
sate does not break translational invariance. Further- 
more, we will show that the mean-field description of a 
well-known model from condensed matter physics which 
is used in the description of spin-Peierls systems — the 
quasi one-dimensional extended Hubbard model at half- 
filling — is equivalent to the massless generalized Gross- 
Neveu model with two coupling constants. In particular, 
we do not have to impose any restrictions on the sym- 
metry of the condensate and the phase diagrams of both 
models can be identified. 

On the one hand, we aim to supplement the current 
literature by elucidating the mechanism of dimensional 
reduction in the context of selected nonrelativistic field 
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theories. On the other hand, we want to explore the 
fascinating applications of the GN model to condensed 
matter physics. Furthermore, since Hubbard models are 
extensively studied both analytically as well as numer- 
ically, our results might point to applications of those 
theories to the GN model. 

This paper is structured as follows: In Sec. II we give 
a brief overview of the basics of the massless GN model 
with two coupling constants. Section III presents the di- 
mensional reduction of BCS theory onto the GN model 
with discrete chiral symmetry. The reduction holds only 
for homogeneous condensates. Limitations for spatially 
varying condensates as well as the cutoff dependence of 
the phase diagram are discussed. Section IV contains the 
dimensional reduction of the quasi one-dimensional Hub- 
bard model onto the generalized GN model. We give an 
introduction to the model and discuss its particle hole 
symmetry. Exploiting this symmetry, the dimensional 
reduction is formulated. A discussion of the cutoff de- 
pendence and comparison with selected experimental re- 
sults follows. The paper is concluded by a summary and 
conclusions in section V. 



II. MASSLESS GROSS-NEVEU MODEL WITH 
TWO COUPLING CONSTANTS 

In this paper, we consider the massless GN model 
with two coupling constants (genGN), i.e. as defined by 
Eq. ([1} with m = 0. While early work on this system was 
carried out by Klimenko [20[, the full revised phase dia- 
gram and particle content were worked out only recently 
by Boehmer and Thies 21 1. 

We are interested in the phase diagram of this model. 
According to Coleman's theorem, any long-ran ge o rder 
in 1+1 dimensions is destroyed by fluctuations [22|]. As 
shown by Witten (23[, those fluctuations are suppressed 
in the limit N — > oo. The counting of various orders, as 
carried out in different contexts by Dolan and Jackiw [24[ 
as well as 't Hooft (25[, reveals that for TV — > oo a finite 
leading order can only be obtained if g 2 N and G 2 N are 
kept fixed. Subleading orders are suppressed by powers 
of N. 

The suppression of fluctuations in the large- N limit 
admits a semiclassical treatment [H [26J. The Hamilto- 
nian of the genGN model in a mean-field approximation 
becomes 
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2g 2 N 2G 2 N 
(2) 

where the scalar and pseudoscalar fields S and P satisfy 
the self-consistency relations 
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-g 2 (^ip) and P = -G 2 ^ 5 ^) 



(3) 



and (. . .) denotes the thermal expectation value. 

The canonical transformation ip — > e^^^tp maps the 
scalar and pseudoscalar interaction terms onto one an- 
other. Hence, we can assume without loss of generality 



that < G 2 < g 2 . Rcnormalization can be performed 
using the conditions [2l[ 



g 2 N 



= In 



v/(A/2) 2 + l + A/2 



In A and 



G 2 N 



= lnA + £, 



(4) 



(5) 



where we set the scalar condensate equal to 1. £ > is 
the renormalized parameter that describes the imbalance 
in the scalar and pseudoscalar coupling. The case £ = 
corresponds to the xGN model whereas for £ — > oo we 
recover the discrete GN model. 

The genGN model contains meson and baryon bound 
states. For a detailed analysis see [27J and [2lJ. The 
phase diagram of the model is shown in Fig. [TJ At zero 
density (small /i), there is a homogeneous phase with van- 
ishing pseudoscalar condensate. A sheet of second order 
critical points separates this phase from a massless homo- 
geneous one. Since in both phases the pseudoscalar con- 
densate vanishes, the critical temperature depends on /i 
but not on £. A sheet of first order transition lines divides 
the massive homogeneous phase from a phase where the 
condensate takes the form of a soliton crystal which inter- 
polates between the kink-antikink condensate of the GN 
model and the helical condensate of the xGN model. All 
transition sheets converge in a line of tricritical points. 
The onset of the inhomogencous phase at zero tempera- 
ture corresponds to twice the baryon mass. Thereby, at 
low density and large £, the condensate takes the form 
of separated kinks and antikinks and goes over into a 
sinusoidal shape for large density. As £ — > 0, the pseu- 
doscalar condensate is less suppressed and the ground 
state oscillates between the scalar and the pseudoscalar 
condensate. It is interesting to note that chiral symme- 
try is never restored at zero temperature. In all cases, 
the spatial period a of the inhomogencous condensate is 
determined by the first gap which opens up at the Fermi 
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FIG. 1: Phase diagram of the genGN model as derived in [2l[ 
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surface: a = n/pp, where pp is the Fermi momentum. 
The fact that a system can lower the energy of its ground 
state by opening up a band gap at the Fermi surface is 
a well-known phenomenon in condensed matter physics 
where it was proposed by Peierls for Id metals 



III. BCS THEORY 

We start with the following observation: The phase di- 
agram of the massless GN model was derived for the first 
time by Wolff Q assuming only translationally invariant 
phases. The diagram is shown in Fig. [2j However, an 
equivalent phase diagram had been obtained by Sarma 
[281 ] even a decade before the proposal of the GN model 
when investigating the phase diagram of BCS theory in 
an external magnetic field [29j . The main difference is 
that the phase diagram of the GN model shows T vs fi 
whereas Sarma derives the critical temperature as a func- 
tion of an external magnetic field. The massive phase of 
the GN model corresponds to the superconducting phase 
of BCS theory and the massless phase with restored chi- 
ral symmetry to the gapless normal phase. This indicates 
strongly that both theories are equivalent on that level. 
This section presents the proof of this equivalence. 

The electrons in a conductor interact via the Coulomb 
interaction with the nuclei and each other. Since the scale 
of the excitation in a conductor is much smaller than the 
binding energies due to Coulomb interaction, the con- 
ductor can be treated as consisting of a gas of weakly 
interacting quasiparticles. We start with the effective 
grand canonical Hamiltonian describing the low-energy 




FIG. 2: Phase diagram of the massless GN model assum- 
ing only translationally invariant phases. The continuous line 
denotes a second order and the dashed line a first order tran- 
sition. The thin continuous lines mark boundaries of metasta- 
bility. 



excitations in a conductor |30j: 

h= y dw ^_z__ M _ MAa 3^_| (v , tV)) 2. (6) 

"0(x) denotes a fermion spinor with components ?/>-(• (x) 
and ipi(x.). m is the effective mass of the quasiparti- 
cles and [i the chemical potential. We introduce an "ax- 
ial" chemical potential \ia which imbalances the two spin 
species: The chemical potential of spin-^ fermions be- 
comes p + [xa and for spin- 4- fermions \x — /M- Physically, 
Ha describes an external magnetic field or the effect of 
impurities in the material. The coupling g > carries di- 
mension EA~ 3 , with E carrying dimension of energy and 
A of momentum. We perform a Fierz transformation to 
rewrite the interaction term: 



i 



(7) 



where C = i<72 ■ The low-energy excitations are described 
by states whose momenta p are close to the Fermi surface. 
This defines a cutoff 



|p| - Pf \ < A/2 < p F 



(8) 



with pf = y/im/i. In this region, we can linearize the 
dispersion relation 



— -m«MIpI 



■Pf)- 



(9) 



The Fermi velocity is defined by vp — Pf /m. All infor- 
mation about the underlying microscopic theory is ab- 
sorbed in vp and A. For this linearization to hold for 
nonzero pa we assume the hierarchy of scales 



HA < UfA/2 < fi- 



(10) 



Renormalization can be performed using a similar con- 
dition as Eq. flU) (U: 



— = ItlvfA, 
9P 



111) 



where we set the energy scale equal to 1 and define the 
density of states (32| 



o 

Pf 
tt 2 vf 



(12) 



In deriving Eq. (fTTj) we apply a large N argument: Since 
we only consider states close to the Fermi surface, the 
fermion propagator G(p) has nonvanishing support in a 
shell around the Fermi surface defined by Eq. ©. This 
means that for |q| ~ 0(pf) the product of two propa- 
gators / dpG(p)G(q + p) is phase-space suppressed by 
a factor of 0(A/pp)- To leading order, this selects all 
"cactus" or "daisy" diagrams. For further reference see 

fum. 
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By virtue of the above argument we analyze J5f in a 
mean-field approximation. Eq. (|6]) then becomes 



H 



d 3 x <^ ^ (v F (\ - iV| - p F ) - fXAO-s) ip 



2 2 g J 



(13) 



where the BCS-condensate A satisfies the self- 
consistency condition 



and (...) denotes the thermal average. 



(14) 



A. Dimensional Reduction 

To keep track of the dimensionality of the various 
quantities we place the system in a box of length L. We 
define the Fourier transform of the spinors by 



1 



5>(p)e ip ' : 



(15) 



Since the Hamiltonian is rotationally invariant, we can 
switch to spherical coordinates and formally replace the 
angular integration by the sum over a number of iV pat 
patches that cover the Fermi surface: 



(2tt)s 



N ¥ 



—y -y 



(16) 



where we set N F = p 2 F L 2 /ix = Anpp / (2ir / L) 2 . Note that 
this is the number of patches that cover the Fermi surface 
if the system is enclosed in a finite box of length L, hence 
iVp a t = N-p. We introduce the notation il>^(p) for a 
spinor with momentum p F +p in the direction of the i-th 
patch. We label the patch opposite to a patch i by —i. 
Eq. (TT3"]) becomes (suppressing the patch labels) 



H 



A, 



pat 



W {VFP- VA<T3)lp+ ^Alp^Clp* 
P ^ 

LA 2 /v F 



-A^ 1 Ctp 



irpg 



(17) 



The factor gp is dimcnsionless. If we perform a partial 
particle-hole conjugation (not changing the spin quantum 
number) for spin- J, particles, 



(18) 



we obtain after normal ordering the Hamiltonian of the 
massless Gross-Neveu model in a chiral basis: 



H 



pat 



E, \ i as, LA 2 /v F 
W 1 \V F p (73 - [I A + A (7i ) V> H 



ngp 



(v F p 7 1 - pa 7° + A) tp + 



LA 2 /v F 
■ngp 



(19) 



where we choose 7 = (7i,7 2 = — i(72 and 7 = (73. This 
establishes the equivalence between both models on a 
mean-field level. Each flavor in the GN model corre- 
sponds to a patch on the Fermi sphere. Spin-^ and -J, 
are mapped onto the right- and left-handed components 
of the relativistic spinor. The coupling gp corresponds to 
(2/ir) g 2 N in the Gross-Neveu model. 

Unlike in a relativistic field theory, the cutoff is a phys- 
ical quantity in BCS theory. Typical cutoffs are of order 
O(10 2 )-O(10 3 ) in units of the mass scale [ll]- In order to 
establish the full equivalence we must assure that for typ- 
ical values of the cutoff the phase diagram of the theory is 
not significantly distorted and we can take A — > 00 with- 
out loss of generality. The Hamiltonian (fH))) is just the 
Hamiltonian of a massive free relativistic Fermi gas with 
single particle energies e(p) = zty^p 2 + A 2 and chemical 
potential pa phis a c-number term. The grand canonical 
potential density is given by (setting v F = 1) 



A/2 



dp 



-A/2 



1 + e 



)} 



A 2 



2g 2 N 



(20) 



Removing the logarithmic divergences by using the renor- 
malization condition Eq. (0 (or Eq. (JTTJ) ) without taking 
the limit A — > 00 we obtain 



(. + .-^5V)}4(hA-i 




A 2 /A 
In , 

2tt 2 



1 A 4 
16tt (A/2) 2 



(21) 



where we subtracted two trivial "would-be" divergences 



1 

8tt" 



A 



2- 



A. 



(22) 



the second one stemming from the infinite fermion den- 
sity of the Dirac sea. The second line of Eq. (|2"Tj) vanishes 
in the limit A — > 00. Fig. [3] shows several phase diagrams 
for cutoff values A = 2.5,4 and 00. As we can see, we 
can take A — > 00 for typical values of the cutoff. 

It is not possible, however, to include inhomogencous 
phases of superconductors (35l . [3o| in the dimensional re- 
duction in the same vein as sketched above. Consider the 
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FIG. 3: Phase diagram of the massless GN model assuming 
only translationally invariant phases for different cutoff values 
(inner to outer) A = 2.5, 4 and oo. The thin continuous lines 
mark boundaries of metastability. 



Hamiltonian Eq. (|13p for a spherically symmetric conden- 
sate, 

l^^AdqDV^dp + qDC^dpl) 
p q 

-iA*(|q|)^(|p-q|)^(|p|), (23) 

where |q| <C |p|. The requirement for dimensional reduc- 
tion is 

|p±q| = |p|±H_p = | P |±|q| (24) 

for all p. This is only possible for q = 0, i.e. a homoge- 
neous condensate. 



IV. HUBBARD MODEL 

In this section we will consider a system from con- 
densed matter physics that displays a one-dimensional 
instability. For a wide class of organic and some in- 
organic materials conduction is essentially restricted to 
one dimension due to their anisotropic structure. They 
can be described as a family of weakly coupled chains. 
This weak transverse coupling allows these materials to 
circumvent the Coleman-Mermin™ Wagner theorem and 
exhibit long-range order. Such phases are characterized 
by a one-dimensional inhomogeneous charge or spin dis- 
tribution and are, therefore, called charge and spin den- 
sity waves (CDW, SDW). These materials are known as 
spin-Peierls systems. For further reference see (37l.[38|. 

As for BCS theory, the phase diagram of the Hubbard 
model does not depend on the chemical potential \x but 



on an external magnetic field h. At low temperature 
and small magnetic field the system possesses a CDWo 
ground state. The CDWo condensate takes the form of 
a plane wave Ae lCJ "' x . The wavevector Q„ is tilted by 
an angle depending on the lattice spacing with respect 
to the preferential direction. At high density, this con- 
densate is modulated in x direction. This is called a 
CDW X phase. The modulation is of order O (Q n )- A 
CDW j, condensate which is modulated in the perpendic- 
ular directions is possible as well, but is excluded for the 
parameter values that we consider in this paper (39j . At 
high temperature a transition to a homogeneous normal 
phase occurs where the gap vanishes. In this section we 
will show that the phase diagram of the Hubbard model 
can be identified with the phase diagram of the gcnGN 
model, Fig. [2] The massive homogeneous phase in the 
gcnGN model corresponds to a CDWo phase, the inho- 
mogeneous phase to a CDW X phase and the massless 
chirally symmetric phase to the normal phase. 

Early work on the relationship between the original 
Hubbard model |4(j with repulsive interaction and the 
chirally invariant Thirring model was carried out by Filev 
[4lj and extended by Melzer [42| as well as Woynarovich 
and Forgacs (43|. Since we are investigating a more 
general model which contains two competing interaction 
terms and admits a nontrivial phase diagram [2l| . there 
is almost no overlap between our work and theirs. 

A. Definition, Symmetries and Anisotropic 
Hopping 

We start with a system of fermions on a hypercubic lat- 
tice with iVi a t sites that are allowed to tunnel (or " hop" ) 
to nearest lattice sites and to interact with their nearest 
neighbors via a spin-dependent repulsive quartic interac- 
tion. The lattice is bipartite, i.e. it can be divided into 
two sublatticcs A and B such that the interaction takes 
place between fermions on different sublatticcs: 

^ = -EE ^(j)vv(i) -/^^>^(j)wcj) 

<y <ij> a j 

- \ E W E V °°' (V4(i)</v(i) - Mi)) 

a,o> ^ <ij> 

x (vd'(j)VV'(j)-n CT ,(j)), (25) 

The sum runs over the spin indices a and lattice sites i. 
< ij > denotes the sum over nearest neighbors and iV M 
the number of nearest neighbors (N^ = 2d for a hyper- 
cubic lattice). Vv(i) is a fermion spinor at site i with 
spin a and n a (j) = (V4(j)Vv(j)) is the mean occupation 
number of a site, tij is the hopping amplitude between 
adjacent lattice sites which we suppose to depend only on 
the lattice direction. As for the BCS theory, h describes 
the effect of impurities or an external magnetic field. The 
most general symmetric form of the coupling is 

V aa > =U C - U s oo' (26) 
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with U c < and < U s < \U C \. A coupling of the form 
—U s (l + a a') is repulsive for spins of the same type and 
favors the formation of an alternating pattern of spin-f 
and -\. fcrmions — a spin density wave (SDW). A cou- 
pling of the form C/ c < enhances an inhomogeneous 
charge distribution — a charge density wave (CDW). 

The Hamiltonian (|25[) is symmetric under a particle 
hole conjugation defined by 



V>! CT (j) it a 
-V>l CT (j) j e b 



V'l.a), (27) 



where Q„ = {■n/a x , n/a y , tt/cl z ) and Xj = 
(Ji a x, 3i a yi3'i a 'z) and au is the lattice spacing in k 
direction. The minus sign ensures that the kinetic term 
is invariant under this transformation for a bipartite lat- 
tice and the factor a in the magnetic term provides the 
invariance of this term. The invariance of the interaction 
term follows from U^ip^(i)ip a (i)U = 1 — ■ip' f _ a (i)ip- cr (i) 
and the symmetry in the summation over a and a' as 
well as < ij >. The momentum decomposition of the 
transformed fields is 



lat j 

VA CT (Q„- P ). 



N, 



(28) 



Note that the transformation reflects the momentum on 
the Fermi surface. For further reference, we note an iden- 
tity that holds for the expectation value of operator bi- 
linears (when 0) [Hj]: 



(V4(p-q)Vv(p)) 



tr{^t( p _ q )^(p)exp[-/3g]} 
tr{exp[-/3#]} 

-(V>! CT (Qn - p) ^- CT (Qn - P + q)>- 

(29) 



For q = the same calculation yields the expression 

(^COVvCj)) + (V'i (r (j)V , -a(j)) = 1- This symmetry im- 
plies half-filling of the system, i.e. half of the number of 
states are occupied: 



£5>£Cj)Vv(j)> = iv lat . 



(30) 



°~ J 



We note that this Hamiltonian is still symmetric under a 
particle hole transformation. The spin-independent part 
of the energy is antisymmetric under a shift by Q n : 



a p a,cr' p.p' ,Q 



xi(p)i(p'-Q)^(p'), 



(31) 



where the single particle energy-spectrum £ a (p) is given 
by 

£ CT (p) = e(p) — ah = — 2ti cosp^d,; — ah. (32) 



e(jp + Qn) = -e(p) 



(33) 



In particular, every point on the Fermi surface of the 
free field theory at half-filling and ft. = is mapped onto 
another point on the Fermi surface under a shift by Q n . 
This is called nesting and Q„ is known as the nesting 
vector. 



B. Dimensional Reduction 

In a quasi one-dimensional system the hopping am- 
plitude in one direction is much larger compared to the 
others: t x ^> t y ,t z . If we assume that the hopping am- 
plitude is proportional to the overlap of atomic orbitals 
this reflects directly the orbital structure and the config- 
uration of the lattice. A slice through the Fermi surface 
is shown in Fig. 21 We linearize the energy around pf- 

e a (p) sa v F (p - Pf) - 2t y cosp y a y - 2t z cosp z a z - ah, 

(34) 

where the Fermi velocity is defined by 



v F 



de(p) 



dp x 



2t x CLx 



(35) 



c — 7r/2a x ,p y — Pz— 



Again we assume the hierarchy h, t yi t z <C VfPf- 

We analyze the Hamiltonian in a mean-field approxi- 
mation 



II 



p,er 



^(p + Q)W(p)A QCT 
(36) 



Q.o" 



The low-energy behavior is determined by the modes 
close to the Fermi surface. For Q only a fluctuation 
term survives in Eq. (|25[) which we will ignore. It is 
customary to neglect scattering into higher states and 
assume Q sa ±Q n . Hence, we obtain the Hamiltonian 

^ = EE-(p)^(p)^(p) + ^E E v^Kp + q) 



Px o 




-7T/S 



Tr/a y 



Py 



FIG. 4: Slice through the Fermi surface of the quasi one- 
dimensional Hubbard model. 
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where we define 

D Q° = E<^(P - Q)V^(P)) and 
p 

a q-' =J2 V ™ ,D <*°'- ( 37 ) 

a' 

Because of the nesting symmetry (|33|) the Fermi sur- 
faces for p x > and p x < can be mapped onto one 
another for small perturbations. In this case, each p x 
value corresponds to two patches on the Fermi surface 
with p x > and p x < 0, respectively. For the right (up- 
per) Fermi surface (i.e. p x > 0) we introduce the notation 

tp^ip) = Vv(l (?;) + (p,0j_)), where 1« points to the i-th 
patch on the right Fermi surface. If we define left-moving 

spinors by ipll(p) = Vv(l W - Qn + (p.Oj.)) the mean- 
field Hamiltonian separates into a sum over all patches. 
The construction is illustrated in Fig. |4j In this notation 



the action of the particle-hole transformation (j2"5)l on the 
spinors reads 

UiPl/rAp)^ = ipl /R _„(-P)- ( 38 ) 

Note that by definition -Dq ct = D-Qa and by the particle- 
hole symmetry ([28]) 

£*Qn+<?o- = D-Q n + q -cr- (39) 

For our choice of couplings U c < and \U C \ > U s > 
it can be shown that the condensate A is only modu- 
lated in x direction j39j |. Such modulations are called 
CDW X phases (as opposed to CDW y phases, where the 
modulation is perpendicular to the conducting direction) . 
We can now split the summation over p into a sum over 
patches on the Fermi sphere and a summation over the 
p^-component. Equation (|36|) becomes 



N -Z^Z^ [VrAp + 1) ^lAp)) I Aq + [-v F p - ah] 6 q . Q { i> La {p) ) 

a p,q x ^ y / N * 

- ^EE V** 1 { D h,+^D Qn+q ^+D Qn+qa D* Qn+qa ,) . (40) 

a, a' q 



We transform the spinors according to 



^Lt(p) 
1pRl(p) 



ipL-rip) 



(41) 



This transformation maps DQ n+qiJ — > aDQ n+qa and the 
c-number term of Eq. (|40[) becomes 



N 

a, a' q 

The condensate becomes 



Q n + q <?'- 



(42) 



AQ„+« t 5(g) - iP(q) and 



(43) 



where we define 



S(q) = u c (D Qn+qt - D Qn+gi ) 

= fE ff ( fl( 3"+f - D -Qn+qv) and (44) 



iP(g) + 

= -y E (• D Q"+9 CT + £) -Q»+?') ■ 



(45) 



The particle hole symmetry implies the relations 

5(g)* = S(-g) (iP(g))* = -iP(-g) (46) 

and we can rewrite Eqs. (|4U|) and (|42p enclosing the sys- 
tem in a box with length L in x direction: 



]V=EEW> + <?) 4» 



. [ '{v F (p + q)-h]S qfi aS(q)-iP(q)\ 
crS(-q) + iP(-q) [-v F P - h] S q<0 J 



E 



L\S(qW 



^lM ) ^ U C LN ' ^ U S LN 
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E 



L\P(q)f 



(47) 



Finally, if we transform the spin- 4. spinors according to 



{MP)J ~* { i>Up) ) 



(48) 



we can recast the theory in the form of the genGN 
model in a chiral basis with coupling constants g 2 N = 
-U c LN/2v F and G 2 N = U s LN/2v F : 

f = E / tetfix) [~h 5 d x -h+ S(x) 7 ° 
+ + 2(-U c LN/2v F ) 
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P 2 {x)/v F 
2{U s LN/2v F y 



(49) 



where we redefine S/P{q) — > L 1 / 2 S/P{q) in order to give 
proper engineering dimension. The additional spin de- 
gree of freedom results in a doubling of flavors in the 
gcnGN model. The transformed condensates are 



S(q) = 



iP(q) 



(J 

X^^p-gJ^GO) and (50) 

a p 

X! ( D Q n +q<? - D -Q n +q<y) 

a 

yEE^-* 5 *))- (51) 



as one would expect from varying the grand canonical 
potential density = — lntr e~^ H / f3L with respect to S 
and P. 

As for BCS theory, the cutoff is a physical quantity. 
The cutoff dependence of the second order transition line 
between massive and massless homogeneous phases has 
already been determined in Sec. III. As pointed out in 
[2~j| , the genGN condensate can be approximated by the 
variational ansatz 

S{x) = 25i cos{2p F x) and P(x) = 2Pi cos(2p F x) (52) 

in the vicinity of the second order transition between in- 
homogeneous and homogeneous massless phase. The sin- 
gle particle energies can be calculated perturbatively us- 
ing almost degenerate perturbation theory which allows 
us to determine the correction to the grand canonical 
potential of the free Fermi gas: 



normal 
normal 

+ MnSf + 2M 12 S 1 P 1 + M 2 2P?- (53) 



Varying this expression with respect to the parameters 
Si , Pi and p p yields two conditions which determine the 
critical values at the phase transition, 



dct M = M11M22 - M\ 2 — and 



where the coefficients Mij are: 



ddctM_ 

dp F 



0. 

(54) 



Mn = PV 



A/2 
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dp 4pp 
2^p 2 



7T 



2tt 





'A 




J 




2 




) 




M22 

M12 = PV 



Mu + - 

IT 



(56) 



A/2 



dp 2pF 
2tt p 2 - p% 



1 



1 



o/3(p-M) 





' A 




\ A 








- In 




) 


2tt V 


2p~F~ ~ _ 









1 _|_ e /3(p+M) 

(57) 



PV denotes a principal value integration. As A — > 00 the 
second lines of Eq. (f55j) equals —\/ir\xi2pF- The second 
line of Eq. ([57]) vanishes in this limit and Eqs. (|55|) - ([57|) 
arc, of course, equivalent to Eq. (129) of [2l|. Examples 
of second order transition sheets are shown in Fig. [5] 
As expected, even for moderate values of A the distor- 
tion of the phase diagram is negligible. We refrain from 
determining the cutoff dependence of the first oder tran- 
sition sheet which would require an extensive numeri- 
cal Hartrce-Fock calculation (4f| H(| and would not yield 
much physical insight. 

The phase diagram of various spin-Peierls com pou nds 
has been measured with high accurac y |3 7l. l47H49l | 
whereas theoretical investigations [HI, |50L |51| have not 
revealed the full phase diagram to the best of our knowl- 
edge. In particular, the first order transition line be- 
tween massive homogeneous (CDWo) and inhomogc- 
neous (CDWz) phase has never been determined. We 
arc now able to exploit our mapping and confront ex- 
perimental data with a full theoretical phase diagram for 
the first time. Fig ure ((6]) shows a fit to data obtained 
by Hase et al. [49[ for the inorganic spin-Peierls system 
CuGeC>3 . The theoretical phase diagram is fit to the 
scale of experimental data. This corresponds to fitting 
the scale parameters Ao and v F which are set equal to 1 
in the analysis of the genGN model. The phase diagram 
of the genGN model was determined in the fj,-T plane 
for fixed values of f = 0,0.1,0.2,0.4,0.8,1.2,2,3,5 and 



4, 5, 00 




2 10 



FIG. 5: Second order phase transition sheet of the genGN 
model for the cutoff values A = 4, 5 and 00. 
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£ = 2 




- 




normal 




- 




CDW 












CDW X 









0.5 


1 1.5 







h [Tesla] 

FIG. 6: Phase diagram of the inorganic spin-Peierls cuprate 
CuGeOa. The fit is done with £ ~ 2. The theoretical curve is 
fit to the scale of the data. The data is taken from [49l |. 



10 and the fit £ = 2 was chosen from this ensemble. 

V. SUMMARY AND CONCLUSIONS 

In this work we explored applications of the Gross- 
Neveu model to nonrelativistic field theories. Starting 
from the striking observation that the phase diagram of 
BCS theory and Gross-Neveu model coincide when re- 
stricting to phases with translational invariance, we were 
able to map BCS theory onto the massless Gross-Neveu 
model with discrete chiral symmetry on a mean-field 
level. We were able to show that the mean-field Hamil- 
tonian of the quasi one-dimensional extended Hubbard 
model, which is widely used in condensed matter physics 
in the description of spin-Peierls systems, is equivalent to 
a generalized Gross-Neveu model with two coupling con- 
stants. In particular, the phase diagrams of both mod- 



els are equivalent including inhomogeneous phases. This 
model was worked out in detail only recently by Bochmcr 
and Thies [2l|. Relying on their results we were able to 
complete the phase diagram of the Hubbard model and 
confront experimental data with the full phase diagram 
for the first time. 

It is interesting to note that although all models have 
been subject to intensive research over the last three 
decades the correspondence of the phase diagrams has 
not been noticed. In particular inhomogeneous phases of 
the Hubbard model had been discussed long before the 
first phase diagram of the Gross-Neveu model was even 
proposed @, HJ. 

This work supplements current efforts to use the ex- 
pertise on quantum field theoretical toy models in the 
study of phases of QCD at high density by providing 
systems that can be dimensionally reduced most clearly 
The equivalence of the Gross-Neveu and Hubbard model 
might lay the ground for further work. In particular, an- 
alytical solutions for phases of the Gross-Neveu model 
with continuous chiral symmetry that do not affect the 
phase diagram might find a physical counterpart fiol . |52| . 
Recently, there has been a lot of interest in baryon scat- 
tering and other dynamical phenomena in the Gross- 
Neveu model [H, [54j • The extensive analytical work over 
the past decades that was devoted to the study of dy- 
namical phenomena in spin-Peierls systems ("sliding" of 
CDWs) [55[ might be useful. 
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